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ABSTRACT 

We investigate the conditions of star formation in the Large Magellanic Cloud. We 
have conducted a survey for water maser emission arising from massive young stellar 
objects in the 30 Doradus region (N 157) and several other Hll regions in the Large 
Magellanic Cloud (N 105 A, N 113 and N 160A). We have identified a new maser source 
in 30 Dor at the systemic velocity of the LMC. We have obtained 3—4 /im spectra, with 
the ESO Very Large Telescope, of two candidate young stellar objects. N105AIRS1 
shows H recombination line emission and its Spectral Energy Distribution (SED) and 
mid-infrared colours are consistent with a massive young star ionising the molecular 
cloud. N 157BIRS1 is identified as an embedded young object, based on its SED and a 
tentative detection of water ice. The data on these four H n regions are combined with 
mid-infrared archival images from the Spitzer Space Telescope to study the location 
and nature of the embedded massive young stellar objects and signatures of stellar 
feedback. Our analysis of 30 Dor, N 113 and N 160A confirms the picture that the 
feedback from the massive O and B-type stars, which creates the Hn regions, also 
triggers further star formation on the interfaces of the ionised gas and the surrounding 
molecular cloud. Although in the dense cloud N 105A star formation seems to occur 
without evidence of massive star feedback, the general conditions in the LMC seem 
favourable for sequential star formation as a result of feedback. In an appendix we 
present water maser observations of the galactic red giants R Doradus and WHydrae. 

Key words: stars: formation - H n regions - Magellanic Clouds - stars: pre-main- 
sequence - infrared: stars - masers 



1 INTRODUCTION 

One of the great unknowns in our understanding of star 
formation concerns the role played by galactic environmen- 
tal parameters. For instance, molecular cloud processes, e.g. 
cooling and magnetic field diffusion, depend on the presence 
of metals, thus it is unlikely star formation is not affected 
by metallicity. The Magellanic Clouds have metallicity and 
density of the interstellar media (ISM) that are lower than 
in the Milky Way. Thus, they provide a star formation tem- 
plate that is more representative of the formation of stars at 
higher redshifts and offer a step closer towards exploring the 
formation of the first generation of stars at zero metallicity. 

Young pre-main-sequence populations have recently 
been identified in Hll regions in the M agellanic Clouds 
l|Nota et alJ I2006J : ICouliermis et alJ 12006^1 . However, rela- 
tively little work has been done in investigating the ear- 
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lier, more embedded stages of star formation in the Magel- 
lanic Clouds. At the distance of the Large Magellanic Cloud 
(LMC, ~ 50 kpc), one is limited to study only the most mas- 
sive stellar embryos that will form massive O and B stars. 
Infrared (IR) observations are an excellent way to inves- 
tigate the properties of such objects, both because they 
are heavily embedded and thus invisible at shorter wave- 
lengths and also because they are surrounded by dense en- 
velopes of gas, dust and ices that are revealed at these wave- 
lengths Ivan Dishoeckll2004lh Recently, a few massive young 
stellar objects (YSOs) in the LMC have been investigated 
Ivan Loon et alJl2005bt I Jones et alJl2005l) . 

Molecular masers (in particular water, hydroxyl and 
methanol) are very bright emission lines that seem to be 
closely ass ociated with the earli est stages of massive star 
formation ide Buizer et al.ll2005l) when the YSOs are dif- 
ficult to detect even at IR wavelengths. Maser emission is 
found in the vicinity of embedded YSOs and thus they are 
powerful beacons of current star formation in a molecular 
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cloud. Several water maser sources have been discovered in 
the Magellanic Clouds (see below). In this contribution we 
describe a new survey of water masing sources in the LMC 
at 22 GHz with the Parkes Telescope. 

In recent years, the 30 Doradus Nebula (N157) has 
become the paradigm of large scale triggered star forma- 
tion. As the largest and nearest extragalactic Hn region 
in the Local Group it offers a rare i nsight into the spatia l 
and temporal properties of starbursts. lWalborn et all l|l99fll ) 
identified the so-called star- formation fronts, the interfaces 
between the energetic outflows from the central compact 
cluster R 136 (first generation star formation episode) and 
the surrounding molecular clouds, in which there is ev- 
idence of ongoing triggered star formation. Water maser 
emissio n has also been detected at various locations in 
30 Dor JWhiteoak et aT] Il983l : Ivan Loon fc Ziilstral l2001al: 
lLazendic et alJl200jL_ The wealth of bright protostars iden - 
tified in the region l|Rubio et al.ll99ij : lBrandner et alJl200lfl 
hints at the possibility that more masing activity hitherto 
undetected might be present. We describe our mapping sur- 
vey of the inner region of 30 Dor, in which we have identified 
4 distinct maser sources, one of them new. 

Three other Hn regions in the LMC, known to 
show water maser emission (NT13, N 105A and N 160A; 
IWhiteoak et alJll983l: Iwhiteoak fc Gardnerlll986h are also 
observed at 22 GHz, as well as most of the previous non- 
detections for water maser emission in the LMC and SMC 
iScalise fc Bradll982t IWhiteoak et al"1ll983h . Since our ob- 
servations, lLazendic et alJ (120021) also reported 22 GHz wa- 
ter maser emission towards N 159. 

The Hn region N113 llHenizelll956T) shows a complex 
structure of Ha bubbles with a rich molecular gas and 
dust m orphology. Several young clus ters are associated with 
N113 iBica. Claria fc Dottoril ll992D and current star for- 
mation activity is occuring wit hin continuum sources in 
the central area o f the nebula ferooks fc Whiteoakl ll997t 
IWong et alJ I2006D . The Hn region N 105 llHenizd fl95a1 
is a complex of evacuated bubbles and dense molecular 
material with several you ng clusters associated with it 
jAmbrocio-Cruz et alJll998l) : current star formation as indi- 
cated by maser activity seems to concentrate in the denser 
central part of the region, N 105A. N 160A is the brightest 
component in the N 160 H II complex jHenizell956ri. Besides 
mase r emission, both an embedde d protostar (iHen ning et alJ 
1998) and a dense molecular core |B olatto et alJl20 00l ) point 
at ongoing star formation . iHevdari-Malaveri et alJ <|2002T) 
and lNakaiima et alJ i2005l) describe the gas morphology and 
stellar content of N 160. 

In order to investigate the relationship between the 
masers and the gas kinematics and star formation activity, 
we have obtained narrow-band Ha images at the AAT of 
N 113, and L'-band images and 3—4 /im spectroscopy at the 
ESO-VLT in N 157B and N 105A. These new observations, 
combined with archival IRAC/Spitzer images are used to 
relate the embedded population to the local gas and dust 
morphology. 



2 OBSERVATIONS 

The new observations used in this paper are described in 
this section and a summary is given in Tabled 



2.1 The 22 GHz survey at Parkes 

The 64 m radio telescope (effective diameter of 45 m) at 
Parkes, Australia, was used from June 30 to July 9 2001, 
with the K-band (1.4 cm) receiver plus autocorrelator back- 
end to observe the 6ib — * 523 rotational transition of ortho- 
H 2 at a rest frequency of 22.23507985 GHz. The observa- 
tions were performed at ~ 22.216 GHz yielding a velocity 
coverage of ~ 860 kms -1 with 0.42 kms -1 channel - . Using 
the dual circular feed, spectra were obtained simultaneously 
in left and right circular polarisation; these were then aver- 
aged. 

The system temperature varied between 120 and 140 K. 
The conversion factor from antenna temperature to flux 
density was 6 JyK -1 . Observing conditions were sometimes 
rather unstable due to clouds. The on-source integration 
time was 20min per pointing. The pointing and focus were 
checked regularly by observing the bright maser sources 
R Doradus and WHydrae (see Appendix^}. The absolute 
flux calibration is accurate to ~20 per cent. 

The obtained spectra were corrected for two distinct 
baseline effects: a low-frequency feature (easily removed with 
a 3 rd -degree polynomial) and an interference signal with a 
frequency of ~ 2.85 MHz. The shape of this interference was 
re-constructed by median averaging the cycles within each 
spectrum, and was then subtracted from the spectrum. 

The observations in the central area of 30 Doradus 
(LHA 120-N 157A) were performed on a double grid system 
(x,y): the primary grid has a separation of one beam, i.e. 
1.3'; the secondary grid positions are obtained by shifting 
the primary grid positions by 39" in right ascension and 
declination. We performed 58 pointings at 47 distinct posi- 
tions around R 136 (some positions deemed interesting were 
observed twice), covering a 5' x 5' area. The pointings are 
shown in Fig.0 superimposed on a Midcourse Space exper- 
iment (MSX) 8 fim image. An extra pointing was performed 
on N 157B, but note that it is offset by 1' from the peak 
emission in the mid-IR. 

The observations of N113 were performed in a 5- 
position dither pattern with 39" shifts. Single pointings were 
performed towards a number of other Hn regions; maser 
sources were detected towards N 105A and N 160A, while 
for 11 other regions no masing source was detected (see Ap- 
pendixIXl . 

2.2 The 3 4 /im observations at the VLT 

The Infrared Spectrometer And Array Camera (ISAAC) 
on the European Southern Observatory (ESO) Very Large 
Telescope (VLT), Chile, was used on 7 and 8 Decem- 
ber 2003 to obtain long-slit spectra between 2.85 and 
4.15 urn of two candidate emb edded YSOs. MSX LMC 888 
fegan. van Dvk fc Pricell200ll) (henceforth N157BIRS1) is 
located in the 30 Dor region and was selected based on its 
very red colours (consistent with a massive YSO) while 
MSX LMC 8 (N105AIRS1) is a candidat e protostar first 
identified bv lEpchtein. Braz fc Sevrel (Il984h — see Sect.EHl 
for a complete discussion on these sources. 

The resolving power of A/AA ~ 700 was set by the 
~ 0.5" seeing rather than the 2" slit width. The thermal-IR 
technique of chopping and nodding was used to remove the 
background, with a throw of 10", jittering within 2" to cor- 
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Table 1. Overview of the new data on the Hll regions discussed in detail in this paper. 
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Figure 1. 22 GHz pointings in 30 Dor superimposed on an MSX image at 8 ^m. On the left, the + symbols identify R136 and the 
published positions of the maser sources 0539—691 and 0539—691 B; o symbols show the positions of the maser sources 0539—691 (a 
and 5), 0539— 691B and 0539— 691C determined in this paper (Section l3.1,ll , On the right, the positions observed are schematically 
represented and numbered. An extra pointing was performed towards N 157B, in the lower right corner of the image. The greyscale is 
logarithmic between 10 — 10 and 3 X 10 — 9 Wcm sr . 



rect for bad pixels. Total exposure times were 12 minutes. 
The spectra were extracted using an optimal extraction al- 
gorithm. An internal Xe+Ar lamp was used for wavelength 
calibration. The relative spectral response was calibrated by 
dividing by the spectrum of the B-type star HIP 020020, ob- 
served within an hour from the targets at the same airmass 
and similar seeing. This removed most of the telluric absorp- 
tion lines — except the saturated methane line at 3.32 /jm — 
but introduced artificial emission features due to the photo- 
spheric lines of HIP 020020. The spectra were therefore mul- 
tiplied by a hot blackbody continuum with Gaussian-shaped 
absorption lines of Bra 4.052, Pfy 3.741, PfS 3.297, Pfe 3.039 
and PfC 2.873 /im scale d to match those i n HIP 020020. More 
details can be found in Ivan Loon et al.l (2006). 

The acquisition images were obtained through the L' 
filter (3.78 /im, bandwidth AA = 0.58 /im). At a scale 
of 0.148" pixel -1 , the (array- windowed) field-of-view was 
57" x 57". The images were obtained in chopping-only mode, 
with a throw of 10" in the North-South direction. The total 
exposure time was 20 seconds. The photometric accuracy is 
of the order of 0.1— 0.2 mag, found from cross-correlation of 
acquisition images of other targets in the same observing ru n 
with dedicated L'-band photometry llvan Loon et all2005aT) . 



2.3 The Ha mapping at the Anglo-Australian 
Telescope 



The Taurus Tunable Filter instrument on the Anglo- 
Australian Telescope (A AT), Australia, was used on 17 July 
2001 to obtain narrow-band Ha images of N113. The blue 
etalon was used in combination with the EEV CCD, at 
a spectral resolution of AA ~7A around Ha. The central 
wavelength was set to Ao = 6568 A, corresponding to a 
Doppler shift similar to that of the LMC. The technique of 
charge-shuffling was applied to take a series of three consec- 
utive images before reading out the CCD, where each image 
was taken on the same area of the CCD but at a central 
wavelength of Ao + AA, Ao and Ao — AA, respectively. Each 
image received an integration time of 1 minute, and a total 
of six such cycles were performed before CCD read-out. This 
was repeated with images taken at Ao + 3AA/2, Ao + AA/2 
and Ao — AA/2, respectively, to improve the spatial and spec- 
tral sampling of Ha; as the instrument is placed in the pupil 
the wavelength changes with radial distance r in pixels as 
(cf. iBland fc Tulivlll989l) : 
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(1) 



where p — 13.5 /im is the pixel size and / = 127.8 mm 
the camera focal length. The two sets of observations were 
repeated with an offset of 1' in the North-South direction 
to cover more of N113. The observations were made at a 
distance of 50° from the zenith and a seeing of 1.8". 

The data were corrected for bias offset and divided by a 
Hatfield image. For each of the two pointings the images were 
combined to produce one image corresponding to the Ha line 
emission and another image corresponding to the underlying 
continuum emission. This was done on a pixel-by-pixel basis 
by assuming that the Ha emission has a gaussian spectral 
shape with a FWHM of 7 A, and that the continuum is 
flat. The resulting Ha image is very clean with most stars 
removed. The main Ha emission structures are so strong 
that they are still present in the "continuum" image, albeit 
at a much lower level fSect. l3.4jl . 



2.4 Additional infrared data 

Near-IR JHK S phot ometry is taken fr om the 2-Micron All- 
Sky Survey f2MASS. ICutri et al.ll2003Tl . Embedded massive 
stars have very red J—K s colours, thus, at the distance of the 
Magellanic Clouds, most embedded YSOs are only detected 
by 2MASS in the K s -band. Hence, bright K s -band sources 
with faint or undetected shorter wavelength counterparts in 
2MASS are embedded YSO candidates. 

For the two luminous IR objects, N157BIRS1 and 
N105AIRS1 (Sect-EOJ, mid-IR photometry at 8.28, 12.1, 
14.7 and 21.3 /im is taken from version 2.3 of the Mid- 
course Space exp eriment (MSX) Point Source Catalogue 
(Egan "et alJ 120031) . The spatial resolution of these data 
varies from 7" at 8.28 ^m, to 18" at 21.3 /im. We also 
collected scans from the IRAS data server 1 for these two 
sources, to measure their 12, 25, 60 and 100 (im flux densi- 
ties using the Groningen Image Processing SYstem (GIPSY) 
software with the SCANAID tool. These measurements were 
fully consistent with the IRAS Point Source Catalogue 
feeichman et alJll988l) . but contrary to the PSC we were 
able to obtain a measurement for N 157B IRS1 also at 60 and 
100 /jm. In addition, for N 105A IRS1 we inspected Low Res- 
olution Spectrograph scans from the IRAS data server. The 
three complete scans are fully consistent with each other, 
and show a weak depression around 10 fim which could be 
indicative of absorption by silicate dust. The mid-IR pho- 
tometry helps constrain the luminosity and hence the mass 
of the star. 

In our analysis of the morphology of the 4 H II re- 
gions we also made use of archival images obta ined with 
the Spitzer Space Telescope JWernere^alJ I2004D . using its 
Infrared Array Camera (IRAC. lFazio et alJ200l . IRAC im- 
ages of 30 Dor were obtained under Early Release Observa- 
tion programme #1032 (P.I. B. Brandl). Both N 105A and 
N 113 observations are part of the cycle 2 Legacy programme 
#20203 (P.I. M. Meixner) while the N 160A observations 
are part of the Guaranteed Time Observations programme 
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Figure 2. Intensity map of the 22 GHz survey of 30 Dor. In this 
contour plot, the darker regions indicate possible masing sources. 
Open circles indicate the grid positions that were observed, cor- 
responding to the pointings in Fig.0 Position (1,1) corresponds 
to (RA,Dec) = (05 h 39 m 11.6 s , -69°09'57"), and the grid spacing 
is half a beam (39"). The intensity scale in this plot is from Jsl to 
5 Jykms - . Comparison with Fig.fflsuggests that maser emission 
seems to concentrate towards the brighter areas of the nebula at 
IR wavelengths. 



#124 (P.I. R. Gehrz). Extended pipeline products (i.e. flux- 
calibrated image mosaics) for the 4 IRAC bands (3.6, 4.5, 
5.8 and 8.0 /im) were retrieved from the Spitzer archive 2 . We 
have performed aperture photometry on the IRAC mosaics, 
for the sources present in the L' acquisition images. We per- 
formed aperture corrections but not ar ray-location depen - 
dent corrections or colour corrections jReach et al] 120051. 
so fluxes might be uncertain by as much as 10 per cent, still 
fully adequate for our analysis (see Sect. l3.3L 



3 RESULTS 

3.1 22 GHz survey 

3.1.1 Mapping of 30 Dor (N 1 57) 

The goal of the 30 Dor survey was to detect maser sources 
and locate their positions to better than the beam size. We 
assume that the sensitivity of the telescope beam is a func- 
tion of the distance to beam centre, represented by a gaus- 
sian distribution. We reconstruct the spectra at each posi- 
tion in the observed grid by combining each spectrum with 
the available adjacent spectra in the grid, weighted by their 



1 http:/ /www. astro. rug.nl/IRAS-Server/ 



2 http:/ /ssc. spitzer. caltech.edu/archanaly/ 
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Figure 3. Maser sources detected in 30Doradus. The spectra cor- 
respond to the positions 0539-691, 0539-691B and 0539-691C 
respectively. The masing structure at 0539—691 (pointing 11) has 
two distinct components (a and b as indicated) at velocities con- 
sistent with the local ISM velocity. The maser at 0539— 691B 
(pointing 3 6) has a strong ly supersonic velocity as observed by 
Ivan Loon fc; Ziilstral l2001ah . The new maser source 0539-691C 
(pointing 25) has a single component at the systemic velocity. 



variance and the gaussian weight corresponding to the dis- 
tance between the grid points. With this process we improve 
signal-to-noise ratio and detection sensitivity. By integrating 
the flux of the strongest spectral feature in each spectrum, 
we construct the intensity map shown in Fig.|2| This map 
shows the location of the strongest intensity peaks; in a sec- 
ond step the spectra at each of those intensity peaks and 
closely surrounding ones are checked (e.g., for velocity con- 
sistency) to separate bona-fide maser sources from spurious 
peaks. The detection sensitivity is ^Uykms" 1 (~5cr). It 
can been seen from this image that there are 2 main flux en- 
hancement regions: a complex region towards the top, that 
actually includes multiple masing sources and a single source 
area towards the bottom of the diagram. 

On kinematic and flux intensity grounds, we have iden- 
tified 4 distinct water maser peaks in 30 Doradus. In each 
spectrum where the maser emission is present, we measure 
the peak intensity of each component at the velocities listed 
in Tabled We obtain a spatial representation of the inten- 



sity of each component that we fit with a two-dimensional 
gaussian distribution with three free parameters: the ampli- 
tude and position (x,y) of the centre — the width of the 
distribution is fixed by the beam size. The amplitude trans- 
lates into the intrinsic intensity of the source and the centre 
is converted into right ascension and declination coordinates 
of each component, with a typical positional accuracy of a 
quarter of a beam size, ~20". In Table|5|we list the mea- 
sured velocity, intensity and position of each maser source, 
as well as their identifications (see below). The last column 
provides the pointings at which the masers were detected 
(Fig-0- The maser spectra are presented in Fig|3] 

0539-691 was first detected at 22 Ghz by 
IWhiteoa ke^tlJ l^83j) , in a regio n ~ 1.5' northeast 
of R 136. Ivan Loon fc Ziilstral fcOOlal) further analysed this 
region, detecting 0539—691 as well as a new water maser 
source, 0539— 691B, located ~ 2.5' southeast of R 136. 
On kinematic grounds the authors distinguish two masing 
sources in 0539—691, one with a velocity consistent with 
the systemic velocity of the local ISM (vh ei ~ 270 kms" 1 ) 
— thi s is the source originall y identified bv | W hiteoa k et alJ 
l|l983l) and also observed bv lLazendic et al.l i2002l) — and 
a component blue-shifted by ~ 70 kms" 1 . 0539— 691B was 
also found to be blueshifted with respect to the systemic 
velocity by ~ 90 km s" 1 . 

We find that the maser 0539—691 comprises 2 compo- 
nents (a and b) with peak velocities ~ 261 and 270 kms" 1 . 
They are very close (~ 1.0') and component b is rather 
wea k, thus it is con ceivable that it is a single source. 



lLazendic et al.l i2002D also detect a component with peak 
velocity 269.5 kms -1 and they pinpoint its location very ac- 
curately. They claim that their spectrum also contains ad- 
ditional faint emission extending down to 258 kms" that 
might be related to our component b. We did not de- 
tect the bl ueshifted compo nent detected in 0539—691 by 
Ivan Loon fc Ziilstral <l2001al) . We do identify a weak source 
with peak velocity of ~ 194 kms" 1 , located ~ 4 . 6' from 
0539-691, with 0539-691B llvan Loon fc Ziilstral lioOlU) . 
We have discovered a new maser source closer to R 136 
(~ 2.8' South from 0539-691), that we call 0539-691C. It 
has a velocity of ~ 266 kms -1 , consistent with the systemic 
velocity. Despite the large survey area, this is the only new 
component we identified. No maser emission was detected 
towards N 157B, at an rms level of 0.13 Jy but note that the 
target was not well centered in the telescope beam. 

3.1.2 Sparse map of N 113 

W ater maser emission at 22 GH z was first observed in N 113 
bv lWhiteoak fc Gardner! dl986Tl . We performed 5 pointings 
in this region. The method used to analyse the spectra of 
N 113 is similar to what was outlined in the previous section 
for 30 Doradus, except with fewer pointings. In Fig.^l we 
plot the central spectrum of our mosaic of pointings in N 113, 
with 4 components identified as a, b, c and d. 

We were not able to disentangle the positions of the 
different components, so we cannot establish positively how 
many spatially-distinct masers are there in this region. How- 
ever, from the analysis of the peak intensity variations 
in our mosaic of pointings, we find evidence that com- 
ponent c might originate from a different position than 
the remaining components: at the position northwest of 
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Table 2. Observed properties of water maser detections: peak velocity, peak intensity, source position and pointings identification (for 
30 Dor only, Fig.0. We estimate that the positions in 30 Dor are accurate to about a quarter of a beam size, approximately 20". 
0539— 691C is a hitherto unknown component. For the other sources we are not able to compute accurate positions; for these, telescope 
pointing positions are listed. No maser emission was detected towards N 157B and thus we provide the rms level in the spectrum. 



Maser source 


peak velocity 


peak intensity 


RA 2000 


Dec 2000 


pointings 




(km s — 1 ) 


(Jy) 


m s J 


(° ' ") 


identification 


0539-691a 


269.5 


2.75±0.60 


05 38 49.9 


-69 04 34t 


6,7, 10,11,12,14, 15,18, 19 


0539-6916 


261.0 


0.38±0.09 


05 38 55.2 


-69 04 12 


6, 10, 11, 18, 19 


0539-691B 


194.0 


0.36±0.09 


05 39 04.3 


-69 07 52 


32,36 


0539-691C 


266.0 


1.01±0.27 


05 38 47.1 


-69 06 06 


19, 20, 25, 29, 30 


N113a 


253.0 


73.8±16.0 


05 13 23.1 


-69 22 34 




N113b 


254.5 


19.4±4.4 








N113c 


251.0 


4.6±1.1 








N113d 


258.0 


2.0±0.5 








N160A 


253.0 


3.3±0.7 


05 39 42.7 


-69 38 26 




N105A 


260.0 


1.5±0.3 


05 09 50.7 


-68 53 23 




N157B 




0.13 (rms) 


05 37 40.2 


-69 11 00f 





f a more accurate position for this source can be found in Lazc ndic et alj l2002fl . 
\ this position is 1' to the West of the peak IR emission. 
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Figure 4. Water maser sources detected in N 113. This spectrum 
is the central one from our 5-position sparce map (see Section f2.ll . 
The different components are labelled a, b, c and d. 



the central pointing, the peak intensity of component c 
was reduced by ~ 60 per cent, while the peak intensi- 
ties of the other componen ts were reduced to ~ 96 per 
cent. lLazendic et alJ l|2002l) recently identified 2 distinct 
masers sources in N 113, 26" apart, and with peak velocities 
253 and 250 km s . We id entify our mult i ple co mponent 
source a, b and d with the lLazendic et alJ 11200211 compo- 
nent at 253 km s^ 1 ; base d on the peak velocity , we believe 
our component c is the lLazendic et alJ <|2002T) component 
at 250 km s " 1 . The stronge s t mas ing component is weaker 
than in the lLazendic et all i2002fl measurements. We note 



that water maser emission associated with YS Os is known 
to be variable both in intensity and velocity iPalaei et alJ 
119931: iTofani et alJll995fl . 

3.1.8 Other regions in the LMC 

Fig.[5]shows the maser detections towards N 160A (top) and 
N 105 A (bottom). We are unable to derive accurate posi- 
tion information from these obser vations but both thes e 
sources have been observed also bv lLazendic et all {2002), 
providing positions with subarcsecond accuracy. The maser 
spectrum of N 105A is similar to that reported previously, 
both in kinematic complexity and intensity. For N 160A two 
maser locations, approximately 44" apart have been iden- 
tified previously at essentially the same velocity. It is very 
likely that both those components contribute to the signal 
in our spectrum; indeed the flux ratios of the brighter com- 
ponent at 253 km s -1 to the other two components at 248 
and 259 kms^ 1 are consis tent with what would be expected 
if the two components of lLazendic et ail i2002T) were com- 
bined. 

3.2 VLT 3 4 fj,m observations 

The L'-band acquisition images of N 157B and N 105A allow 
us to identify a number of red sources other than the two 
spectroscopic targets. Aperture photometry was performed 
on all sources in each image. Very red objects (typically 
K s — L' > 1 mag) are referred to as "IRS" , and these are 
numbered on the basis of the K s — L' colour and L'-band 
brightness. These IRS objects have significant excess emis- 
sion at 3 /im (Fig.[7J and are not simply reddened objects 
behind the molecular clouds. The other point sources are 
referred to as "S", and these are numbered with increasing 
RA. Two other sources in N 105A seem to be diffuse in na- 
ture, at least in the L'-band, so they are named blob A and 
B (see discussion below). The identification, 2MASS posi- 
tions and magnitudes and L'-band magnitudes are listed on 
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Table 3. List of objects in N 157B and N 105A for which we have 3— Afim data, in order of increasing RA (all J2000 coordinates are 
based on 2MASS). Near-IR magnitudes are from 2MASS (JHK S ) and our L'-band acquisition images. 



T T T A 1 in 

LHA 120-.. . 


HA 2UUU 


Dec 2000 


J 


H 




L 


J— K a 


K s — L 




fh m 


(°, ' arcsec) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


N157BS1 


05 37 45.5 


-69 11 10 


13.14 ±0.05 


12.07 ±0.04 


11.58 ±0.03 


11.11 ±0.01 


1.56±0.06 


0.47±0.03 


N157BS2 


05 37 46.0 


-69 11 15 


> 13.2 


> 12.6 


12.95 ±0.05 


12.45 ±0.06 


> 0.3 


0.50±0.08 


N 157B IRS2 


05 37 46.1 


-69 10 50 


> 17.0 


> 15.6 


14.33 ±0.09 


11.39 ±0.02 


> 2.7 


2.94±0.09 


N157BS3 


05 37 46.4 


-69 11 17 


13.35 ±0.04 


13.03 ±0.06 


12.89 ±0.06 


12.76 ±0.07 


0.46±0.07 


0.13±0.09 


N 157B IRS1 


05 37 50.3 


-69 11 07 


15.93 ±0.12 


14.22 ±0.07 


11.45 ±0.03 


7.57 ±0.01 


4.48±0.12 


3.88±0.03 


N 157B S4 


05 37 51.4 


-69 10 56 


13.81 ±0.05 


12.79 ±0.03 


12.50 ±0.04 


12.27 ±0.05 


1.31±0.06 


0.23±0.06 


LHA 120-.. . 


RA 2000 


Dec 2000 


J 


H 


K s 


L' 


J-K s 


K s -L' 




m s~j 


m s'j 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


N105AS1 


05 09 45.6 


-68 53 19 


14.31 ±0.03 


13.41 ±0.03 


13.28 ±0.04 


13.31 ±0.16 


1.03±0.05 


-0.03±0.16 


N 105A S2 


05 09 47.5 


-68 53 36 


11.68 ±0.02 


10.87 ±0.02 


10.63 ±0.02 


10.41 ±0.01 


1.05±0.03 


0.22±0.02 


N 105A S3 


05 09 48.4 


-68 52 51 


14.65 ±0.08 


13.84 ±0.04 


13.70 ±0.09 


13.43 ±0.14 


0.95±0.12 


0.27±0.17 


N 105A S4 


05 09 48.5 


-68 52 44 


13.89 ±0.05 


13.08 ±0.04 


12.98 ±0.05 


12.78 ±0.07 


0.91±0.07 


0.20±0.09 


N 105A S5 


05 09 49.1 


-68 52 54 


14.67 ±0.05 


13.90 ±0.07 


13.89 ±0.07 


13.74 ±0.17 


0.78±0.09 


0.15±0.18 


N 105A S6 


05 09 50.3 


-68 53 35 


15.42 ±0.09 


14.06 ±0.06 


13.73 ±0.07 


13.17 ±0.12 


1.69±0.11 


0.56±0.14 


N105AIRS1 


05 09 50.6 


-68 53 05 


> 15.3 


> 14.7 


13.77 ±0.10 


9.88 ±0.01 


> 1.5 


3.89±0.10 


N105AIRS2 


05 09 52.0 


-68 53 29 








14.30 ±0.34 






N105AS7 


05 09 52.3 


-68 53 20 


> 15.9 


15.25 ±0.14 


14.16 ±0.15 


13.61 ±0.19 


> 1.7 


0.55±0.24 


N 105A blob A 


05 09 52.9 


-68 53 01 


14.53 ±0.09 


13.88 ±0.11 


13.02 ±0.10 


12.83 ±0.11 


1.51±0.13 


0.19±0.15 


N 105A blob B 


05 09 52.5 


-68 53 02 














N 105A S8 


05 09 53.8 


-68 52 53 


13.40 ± 0.04 


13.21 ±0.04 


12.88 ±0.05 


12.35 ±0.10 


0.52±0.06 


0.53±0.11 


N 105A S9 


05 09 55.6 


-68 53 13 


14.39 ±0.04 


13.34 ±0.04 


13.11 ±0.03 


12.96 ±0.17 


1.28±0.05 


0.15±0.17 



Tabled In the next subsections these sources and the IR 
spectra are described in detail. 

3.2.1 Infrared sources in N 157B 

MSXLMC 888 (=IRAS 05381-6912. lEgan et aljElOlh is a 
bright mid-IR point source in the 30 Dor region, situated 
within 1' fromN157B (usu ally identified as SNR0538-69.1, 
although lChu et alJ l)2004Fl question the nature of the nebula 
as a supernova remnant). Close to N 157B there is a small 
molecular cloud, JGB30Dor-22, w ith a diameter of 6.8 p c 
identified from CO observations l|,Tohansson et all fl998'l. 
The cold emission from the molecular cloud appears spa- 
tially resolved on MSX images, whilst the bright mid-IR 
source stands out through its warmer unresolved emission. 

In our L'-band image (FigJfjJ, the MSX source is eas- 
ily identified as a bright IR star, that we call N157BIRS1. 
There is at least one other star nearby also with extreme 
K B - L' colour, N 157B IRS2. IRS2 is seen in the direction 
of the molecular cloud and may be located behind it, but 
its very extreme K s — L' colour suggests it could also be a 
young object embedded in the cloud, warming the surround- 
ing dust. IRS2 is considerably fainter than IRS1. N 157B SI 
and S4 have (J - K s ) > 1.3 mag but (K s - L') < 0.5 mag 
and are likely suffering from extinction by dust within the 
molecular cloud. 

The spectrum of N157BIRS1 (Fig-HJ seems at first 
glance featureless. It does not resemble spectra of evolved 
stars, e.g. a heavily reddened, late-M t ype giant (cf. 
iMatsuura et alJ 120051 : Ivan Loon et alJ [2006) . There is also 
no evidence for H recombination emission lines arising from 
an ionised region (Sect. lpT^^l . When the spectrum is plot- 
ted as optical depth with respect to a pseudo-continuum 
(a low-order polynomial fit) there seems to be a hint of 



the broad water i ce feature at 3.1 /xm, s imilar to that in 
IRAS 05328-6827 ivan Loon et al.ll2005bl) also in the LMC. 
The evidence is, however, not conclusive with a column den- 
sity at most AT(H2 0) < 10 17 cm -2 . It is unclear what are 
the narrow absorption features blueward of 3.5 /im. 

3.2.2 Infrared sources in N 105 A 

MSXLMC 80 f=IRAS 05101-6855. lEgan et alJliool is a 
bright IR source in the Hn region LHA 120-N 105A. It 
is associated with an IR ob ject that was suggested by 
lEpchtein. Braz fc Sevrd <ll984l) to be a "protostar" . This ob- 
ject, which we shall refer to as N105AIRS1, is very bright 
in the L'-band (Fig.|ljJ and with K s — L' = 3.9 mag it is ex- 
tremely red (Table|!|Jl. It is not detected in the J and H-band 
of 2MASS, so it is probably heavily extincted. 

LHA 120-N 105A blob is associated with a small but ex- 
tended nebulosity at a projected distan ce of 14" from IRS 1 
and 10" from S8, the WR star Br 16a fctopita et alJll994l) . 
We identify a bright core, blob A, surrounded by patchy 
emission visible in the L'-band image (Fig©. This diffuse 
core has a point-source 2MASS counterpart, but the 2MASS 
images seem to show some diffuse emission too. We have as- 
signed an identifier to one of these patches, blob B, because it 
happened to fall in the slit of the spectrograph and produce a 
clear sig nal (see below). Contin uum emission was detected at 
6.6 GHz fellingsen et al.ll994l) . centred on IRS1 and blob A, 
connecting the two sources. It thus appears that both these 
sources are embedded in an ionised environment (see below) , 
probably created by IRS1 itself or a source within blob A. 
The (K s — L') colour of the bright knot blob A is not partic- 
ularly red (Fig.|7J even though (J — K s ) = 1.51 mag suggests 
some extinction. 

IRS 2 is a very faint source located within ~ 0.5" from 
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Figure 5. Water maser sources detected towards N 160A (top) 
and N105A (bottom). 



the posi tions of the water jLazendic et alJl2002t) and OH 
masers ferooks fc WhiteoakHl997l) . It was not detected in 
the 2MASS survey, but appears as a very faint star in our 
L'-band image (Fig. [(J. Given the low density of stars in the 
image the association with the masers is suggestive. This 
object co uld be an embedded Y SO. The methanol maser de- 
tected bv lSinclair et all lll992l) unfortunately falls just out- 
side the SE corner of the field of the L'-band image (see also 
lEllinesen et~aHll994) . 

Several more sources can be identified from Fig.|S| They 
all appear stellar in nature and have (K s — L') < 0.6 mag. 
S6 and S9 have (J — K s ) > 1.2 mag thus they could be sit- 
uated behind the molecular cloud. The remaining sources 
have unremarkable near-IR colours. 



IRS 2 



S4 



IRS 1 



S 1 




IRS 2 



S 2 



S 6 



Figure 6. L'-band acquisition images of N 157B (top) and 
N 105A (bottom), with identified stellar sources ("S") and in- 
frared sources ("IRS"). For N157B, the slit was centred on IRS1 
and aligned N— S, while for N 105A the slit orientation was chosen 
to include the diffuse source blob A (see text). The faintest stars 
in the images are L = 12.76 mag (N157B S3) and L = 14.30 mag 
(N105AIRS2), respectively. 



3.2.3 Hydrogen line emission and outflows in N 105 A 

The spectrum of N 105AIRS1 shows a very red continuum, 
very likely free-free emission (see Sect . 13.3} . underlying H re- 
combination line emission (Fig.OU. There is no evidence in 
the spectrum for broad ice features. The Bra line is about 
eight times as bright as the continuum, but the Pf<5 line 
cannot be positively detected. This suggests a high dust 
column density in front of the region in which the line 
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Figure 7. Absolute L'-band magnitudes versus K s -L' colours 
for objects near N 157B (open circles) a nd N 105A (filled circles) . 
Approximate empirical sequences from Ivan Loon et al.l l2005afl 
are overplotted for the red giant branch (RGB), red supergiants 
(RSG) and asymptotic giant branch (AGB, with mass-loss evo- 
lution); early type main-sequence stars have K s — L' ~ 0. The 
effect of Ay = 30 mag extinction is also shown. Only three of the 
sources in this diagram are particularly red. 
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Figure 9. Spectra of the massive YSO candidate, N105AIRS1, 
and the blobs N105AblobA & B. The strong 3.32 /im telluric 
methane feature has been removed. Besides line emission from 
Bra and Pf/7 and 8, IRS1 has a spatially unresolved (very) red 
continuum whilst blob A & B display broad 3.3 (im emission at- 
tributed to PAHs. 




A [Mm] 

Figure 8. Spectrum of the massive YSO candidate N 157B IRS1 
in the 30 Dor region. The strong 3.32 fim telluric methane feature 
has been removed. The spectrum shows no H recombination emis- 
sion lines. The lower panel shows the optical depth with respect 
to the continuum (dotted line in the top panel); there might be 
a hint of the broad water ice feature at 3.1 /im. 



emission arises. The Bra line emission is blended at this 
resolution with emission from He I at 4.0378, 4.0410 and 
4.049 (im; this emission is ve ry weak when compared to Bra 
iDrew. Bunn fc HoardlrQQot) so it has a negligible effect on 
the measured fluxes. We have measured the line fluxes of 
Bra and Pfry; the Bra to Pfry flux ratio is ~ 9. This ratio 
can be compared to the prediction from recombination the- 
ory. W e consider case B treatment (e.g.. iHummer fc Storevl 
Il987fl and adopt T c = 7500 K and n e = 1000 cm" 3 — the 
emission coefficients a re not very sensitive to t emperature 
or density. Using the IHummer fc Storevl (^987) computed 
values, the theoretical flux ratio is predicted to be ~ 7.5. If 
we adopt a sim ple power law for the extinction as a function 
of wavelength (Mathis 1990), we estimate that the Bra to 
Pf/y flux ratio for IRS1 is consistent with case B scenario 
affected by Ay ~ 40 mag. This agrees very well with the 
reddening needed to fit the object's Spectral Energy Distri- 
bution (SED, Sect .1331. 

The Bra line is marginally resolved in IRS1, both spa- 
tially and kinematically. To study the spatial distribution 
of the Bra emission we plot a position-velocity diagram, in 
which one axis is the dispersion direction and the perpen- 
dicular axis is the spatial direction along the slit (Fig. llOII . 
A spatial profile for the continuum emission is obtained by 
averaging many columns on the array at either side of the 
emission line. After subtracting this continuum (solid lines 
in Fig. 1101 , the line profile has broad wings and shows a clear 
asymmetry with emission in the WSW moving towards us 
and emission in the ENE moving away. This is strong ev- 
idence for a bipolar outflow originating near the centre of 
IRS1, at a projected velocity of ^bipolar ~ 100 — 200 kms" 
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Figure 10. Position-velocity diagram of the massive YSO, 
N105AIRS1, showing emission in the Bra line before (dotted) 
and after (solid) subtracting the continuum emission. The con- 
tour levels range from 6 x 10 -6 to 1.9 x 10 _4 W m 2 /im 1 sr 1 
in increments of a factor 2. The blueward wing of the Bra line 
profile is affected by blends with very weak He I emission, ap- 
proximately at —1070, 835 and — 245 kms - , responsible for the 
emission tail towards large blueshifted velocities. There is a clear 
asymmetry in the Bra profile (highlighted by the dashed line), 
evidence for a bipolar outflow. 




-800 km/s 



relative velocity 



Figure 11. Position- velocity diagram of the N 105A blob A+B 
nebulosity, showing emission in the Bra line. The contour levels 
range from 1.1 X 10~ 5 to 3.5 X 10~ 5 W m -2 /im _1 sr _1 in in- 
crements of 3.5 X 10 _e W m — 2 /im -1 sr — 1 . Blob A resembles an 
expanding, asymmetric shell whereas blobB is a compact, kine- 
matically cold cloud. 



at a few tenths of an arcsecond (~ 0.1 pc) at either side of 
the star. 

Blobs A and B also show emission in the H recombina- 
tion lines but in this case also PfJ is detected fFig.l9fr. in- 
dicating that extinction by dust is less severe than towards 
IRS1. Continuum emission in this source is extremely faint, 
if at all present. Both components show evidence for broad 
emission around 3.3 /im underlying the Pf<$ line. The shape 
of this emission is consistent with the 3.28 /mi unidentified 
IR feature that is usua lly attributed to Polycycli c Aromatic 
Hydrocarbons (PAHs, lAllamandola et al.lfl985l) . The Bra 
to Pfy flux ratio for bl ob A is ~ 5, w hich might suggest that 
Bra is optically thick iDrew et al]ll993h . 

The spectrum of blob A seems broader than IRS1 and 
exhibits a double-peaked line profile (FigJHj, suggestive of 
an outflow also emanating from this nebular object. The 
position-velocity diagram of the Bra line (Fig. lllH shows a 
kinematically broad emission suggesting a compact outflow 
at a projected velocity of u ou tflow ~ ±300 kms -1 . The spec- 
trum shows departures from spherical symmetry: the bright- 
est emission, at ~ 0.2" to the WSW of the core of blob A, is 
associated with gas that is preferentially moving towards us. 
Blob B, on the other hand, is kinematically cold and shows 
no signs of outflow. This seems to suggest that there might 
be a source powering an outflow in blob A, even if there is 
no evidence of free-free continuum from an ionised medium. 



3.3 Spectral Energy Distributions of N 157B IRS1 
and N 105AIRS1 

We performed aperture photometry in the IRAC images for 
the objects identified i n the L'-band imag es, obtaining both 
fluxes and magnitudes jReach et all2005ri . The 3.6 /im mag- 



nitudes agree quite well for most sources in N 105A and 
N 157B, except for the objects in areas where the "nebu- 
losity" appears denser CFigs. lTH and 1151 . In these areas the 
background is not correctly estimated, and measured fluxes 
are depressed. Nevertheless, the colours of these objects 
can still provide another hint on their nature. We briefly 
discuss here only the objects with measurements in the 4 
IRAC bands: N105AIRS1, IRS2, blob A+B and S2, and 
N 1 57BIRS1 and IRS2 

Ijones et alJ <l2005l . and references therein) defines 
the locus for red giants, Class II YSOs (objects with 
disks), Class I YSOs (embedded objects) and Hn regions. 
N 105A S2 has colours consistent with it being a red gi- 
ant. N105AIRS2, N157BIRS2 and N105A blob A+B have 
colours consistent with Hll regions, i.e. regions with ionised 
gas and less dust column density. N105AIRS1 has colours 
consistent with a very reddened compact Hll region, in 
agreement with the presence of H recombination lines in 
its spectrum, reddening determination and the findings of 
the SED analysis (see below). N157BIRS1 has colours con- 
sistent with an embedded (Class I) YSO. 

The spectral energy distributions (SEDs) of 
N105AIRS1 and N157BIRS1 were reproduced 
with the dust radiative transfer model dusty 
Jlvezic. Nenkova fc ElitzuJll99Sl) . as an attempt to derive 
approximate quantitative information about these YSOs. 
The fits are shown along with the observed SEDs in Fie. 1121 
and the main fit parameters are summarised in Tabled In 
both cases we could obtain a reasonable fit to the SED with 
a radial dust density profile p(r) = pin(r /r[ n )~ 15 (steady 
accretion), and a standard MRN grain size distribution 
(Mat hisl Rumpl fc Nordsiecklll977t) . The dust envelope was 
assumed to be spherically symmetric; this is probably not 
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Table 4. Mid-IR photometry (in Jy) of the embedded YSOs N157BIRS1 and N105AIRS1: IRAC 3.6, 4.5, 5.8 and 8.0 fun, MSX 8.28, 
12.1, 14.7 and 21.3 (im, IRAS 12, 25, 60 and 100 fim. IRAC flux densities are uncertain by up to 10 per cent. 



Name 




IRAC 






MSX 






IRAS 






3.6 fim 


4.5 Aim 5.8 /an 


8.0 Aim 


8.28 /zm 


12.1 Aim 14.7 Aim 


21.3 Aim 


12 Aim 


25 Aim 60 Aim 


100 Aim 


N157B IRS1 


0.173 


0.313 0.483 


0.635 


0.96±0.04 


1.62±0.09 1.89±0.12 


4.2±0.3 


4.0±0.3 


25±3 170±20 


140±40 


N105A IRS1 


0.028 


0.063 0.181 


0.473 


1.72±0.07 


4.70±0.24 6.61±0.40 


20.7±1.2 


8.0±0.3 


56±2 310±10 


400±20 



Table 5. Input parameters for, and results from the modelling with the DUSTY code of the spectral energy distributions of the embedded 
YSOs N157BIRS1 and N105AIRS1. Both make use of a standard MRN grain size distribution (see text). 



Name 


central source 


dust type 


Tin 


Pin 






r ou t 


A v 


L 








(K) 


(gem" 


3 ) 


(AU) 


(pc) 


(mag) 


(L ) 


N157B IRS1 


25 000 K blackbody 


silicate 


730 


8 X 10" 


-19 


350 


1.7 


23 


7.8 x 10 4 


N105AIRS1 


free-free (911 A. ..10 cm) 


0.7 silicate - 


- 0.3 carbon 350 


2 x 10" 


-1!) 


2500 


1.2 


48 


2.5 x 10 s 
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Figure 12. Spectral energy distributions of N157BIRS1 (top) 
and N 105A IRS1 (bottom). The IRAS data have poor spatial res- 
olution and probably include emission from cold dust surrounding 
the YSO candidate. Fits obtained with the DUSTY code are over- 
plotted as dotted curves (see text). 



strictly true. The bolometric luminosity was derived from 
scaling the model SED to fit the observed SED and using a 
value for the distance to the LMC of 50kpc. 

The main difference between the two objects is in the 
nature of the central source. In the case of N 157B IRS1 we 
place a black body of 25,000 K in the centre of the dust 
envelope. Such hot star is required to reproduce the SED 
throughout the JHK bands where the star becomes notice- 
able; at longer wavelengths (including the 3—4 Aim band) 
the SED is dominated by dust emission. Based on the SED 
analysis alone, we cannot entirely rule out the possibility 
that this is an evolved dust-enshrouded red supergiant of 
around 10— 15 Mq. However, this would imply a rather long 
(~ 20 Myr) time-gap between the formation of this massive 



star and the star for mation that is currently ta king place in 
the 30 Dor region — IWalborn fc Blades! < 19971) identified 5 
different stellar populations in the 30 Dor nebula, with ages 
in the range < 1 and 10 Myr. A chance coincidence of such 
rare object with N 157B is extremely unlikely. Furthermore 
the cautious identification of water ice (Sect . 13*72.11 argues 
against it. 

In the case of N105AIRS1, the 3— 4 Aim spectrum 
clearly indicates the presence of an ionised region in- 
side of the dust envelope. We thus represented the cen- 
tral source by a free-free emission object. There can be 
little doubt that this is a young object and the sur- 
rounding dust is not produced by the central star, and 
we there fore used a mixture of amorphous oxygen-rich 
silicates "Ossen kopf. Henning fc Mathisl Il992t) and carbon 
'Hannerlll988l) . The SEP* is rather well reproduced, espe- 
cially throughout the 3—4 /im band, and the optical extinc- 
tion used in the model is consistent with that estimated 
from the hydrogen emission line ratios (see Sect. 13.2*"*!") . 
The dust envelope is cold, again consistent with a cavity 
created by the emerging Hn region powered by a mas- 
sive YSO. The bolometric luminosity of N105AIRS1 sug- 
gests a central star of spectral t ype 05 and about 43 Mq 
"Hanson. Howarth fc Contill997|) . 

3.4 Ha map of N 113 

The Ha images of N 113 (Fig.|nj show a wealth of detail. 
The line emission is so strong that it is still apparent in the 
pseudo-continuum image (though at a much fainter level). 
The main extended structures are three regions of Ha emis- 
sion to the east, north and west of the centre of the Hn 
region: N 113F, N 113C and N 113D, respectively. This is em- 
phasised by an obscuring lane running across the H II region, 
onto which are projected (from east to west) N 113A, N 113B 
and N113E: "a chain of three small intense knots of nebu - 
losity which, together, make up NGC 1877" dHenizelll956h . 
There is a fourth such knot — which we name N113G — 
in between N 113A and N 113F. Although fainter, it stands 
out well against the dark dust lane. 

N 113 C and the low er surface brightness nebula 
BSDL945 "Bica et al.lll999l) to the east of N 113F each con- 
tain a particularly striking example of a shell with a central 
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Figure 13. Optical and IR images of N 113. The Ha emission (top left) and pseudo-continuum (top right) images cover approximately 
4.4'x 6.7' centered at ~ 05 /l 13 m 25 s , -69°22'26". The greyscale is logarithmic between 10 4 and 3xl0 5 Rayleigh for the Ha emission 
image and between 200 and 3000 R for the continuum image. The most important gas structures and stars are identified and labelled. 
The 8/im IRAC/Spitzer image is shown at the bottom left (4.7' x 4.7', c entered at ~ 05''13"'27 s , -69°22'28"; the scale indicates 1'; the 
greyscale is li near between —5 a nd lOOMJysr" 1 ). Maser sources (circles: lLa,zendic et al.l2002tlBr"ooks fc Whiteoakll997l) and a molecular 
core (square; IWong et afll200ft are identified. At the bottom right is a false colour image where Ha emission and stars appear in blue 
while red is molecular cloud material bright in the 8 fim image. No attempt was made to remove artifacts from the IR image, for instance 
the ghost impressions caused by the brightest stars. In all images North is to the top and East to the left. 



star (indicated in Fig. ll3l as "shell 1" and "shell 2"). These 
shells are likely to be stellar wind-blown bubbles with a dy- 
namical timescale of only a few ~ 10 3 yr (for a wind speed 
of order 10 3 kms -1 ). The ce ntral star of sh ell 1 can be iden- 
tified with a B0-0.5III star <Wilcotdll994h . but the central 
star of the larger shell 2 remains anonymous. 

Several OB stars in N 113 have been described in the lit- 
erature. Near the edge of the Ha emission, the evolved B2[e] 
supergiant HD 269217 (= Hen S 89, R82, IRAS 05136-6925, 
MSXLMC216) has enjoyed considera b le att ention. It was 
first mentioned in iMerrill fc BurwelH <ll933T) as a B star 
with strong Balmer line emission. Descended from a 30 Mq 
main sequence star, it has been associated with circum- 
stellar dust and found to move at a heliocentric speed of 
v* ~ 240 kms" 1 dZickgraf et alJll98rj) . Another, little stud- 
ied emission-line star HD 269219 (= Hen S 90) is found fur- 



ther north and away from the H II region. Closer to the core 
of the Hn complex, HV2377 is an M-type suspected su- 
pergiant. Both around and embedded within N 113 are a 
number of Q9— B0. 5 stars — either on the main se quence or 
slightly evolved ferunet et al]|l975HWilcotslll994l) . 



4 DISCUSSION 

4.1 Triggered star formation in the mini-starburst 
30 Doradus and N 157B 

The central area of the 30 Dor nebula (N 157A) exhibits 
a complex star formatio n history. As already mentioned, 
IWalb orn fc Blades] 119971) disentangled five distinct stellar 
populations in the central area of 30 Dor. The population 
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Figure 14. IRAC/Spitzer 8.0 (im image of 30Doradus (N157A), 
with N 157B at the bottom right of the image (North to the 
top, East to the left; the scale indicates larcmin). The greyscale 
of the image is linear between —5 and 150MJysr — 1 . Large cir- 
cles are the maser detections presented in this paper (the size of 
the circle represents the approximate positional accuracy) while 
small circles are maser positions fromlLazendic et al,l J2002L wa- 
ter maser) and iBrogan et alj 120041 . OH maser located slightly 
to the south and west of the water maser position) . Tri angles 
are embedded obje cts (see text) from lBrandner et alJ feOO ifl and 
iRubio et alJ il998F) . The large outline represents the area covered 
by our maser survey. The cross indicates the position of R 136 
(05 ,l 38 m 42 s , -69°06'03"). 

located mainly in the molecular filaments to the west and 
northeast of R 136 shows spectroscopic evidence of extreme 
youth (< 1 Myr). These infant objects are spatially related 
with bright, compact IR sources within a complex nebular 
and dust morphology (Hunter et al.lll995l): pillars of molec- 
ular gas and dust and dark globul es jWalborn et~aHll999t 
IWalborn. Maiz-Apellaniz fc Barb3l200^ . which are being 
photo-evaporated by th e intense radiation f rom the massive 
compact cluster R136 faunter et al.lll995l) . Together with 
the maser sources, this constitutes strong evidence that star 
formation is ongoing in the nebula. 

Fig.lTHshows the 8.0 /im IRAC/Spitzer image of 30 Dor. 
At this wavelength, the extended emission is dominated by 
hot dust and PAH emission. This figure also shows the loca- 
tion of young, embedded IR sources ( J — H and H — K 
> 1.5 mag) identified in the literature jRubi o et alJ Il998t 
iBrandner et al.ll200lHMaercker fc Burton! 120051) as well as 
the maser sources. The water masers in 30 Dor are lo- 
cated in similar structures as the IR sources, surrounding 
R 136, suggesting that current star formation (not just re- 
cent star formation) is also concentrated near the interfaces 
of the molecular cloud complex within the influence sphere 
of the previous generation of massive O and B-type stars. 
The large velocity difference between maser source com- 
ponents (Sect. l3~l.H implies location in distinct structures 



in the ISM. As proposed by Ivan Loon fc Ziiistral fcOOlab . 
the highly supersonic velocities of some of the sources and 
their location near the rim of large gas superbubbles (e.g., 
IWang fc~H clfand 199~ ) strongly suggest that the masers oc- 
cur at the collision fronts of rapidly expanding bubbles of 
ionised gas with the surrounding dense neutral material. If 
the masers trace the velocities of protostars within these 
environments, then this constitutes evidence that star for- 
mation was triggered at different locations by the massive 
stars' feedback. 

Our maser observations survey the central part of the 
30 Dor nebula, covering the immediate neighbourhood of 
R136 and the densest regions of the nebula. Even though 
the survey area is much larger than previous observations, 
we only identified one new maser source. We can compare 
our survey result s to the galactic water maser luminosity 
distribution from IValdettaro et alJ a2001T) . Assuming a me- 
dian distance to high-mass star forming regions in nearby 
spiral arms of a few kpc, we would expect the bulk of the 
water masers in 30 Dor to have integrated fluxes of order 
0.1 — 1 Jykms -1 , with a steep decline at the higher end. 
This suggests that, although close, our survey is not yet 
deep enough to reveal the bulk of the water maser popula- 
tion in 30 Dor. Nonetheless, as all detected water masers are 
located at or near compressed interface regions, our survey 
provides evidence that most of the on-going massive star 
formation occurs in the presence of feedback. 

The situation is less clear for the molecular material in 
N157B, still in the 30 Dor region but located ~7' (100 pc) 
from R136. We identify N157BIRS1 as a candidate pro- 
tostar from the analysis of its SED and from the tentative 
detection of water ice. We have detected no maser emission 
towards N157BIRS1. Despite the 0.9' mismatch in posi- 
tion, our experience with the 30 Dor mapping (see for in- 
stance analysis of 0539-691C) indicates that we would have 
detected emission from this source if stronger than ~Uy. 
We cannot exclude the presence of weaker emission, how- 
ever water masers are normally associated with the earlier 
(hot core) stage when such objects are very weak in the IR 
ide Buizer et al.ll2005l) — IRS1 is the brightest source in 
our sample at near and- mid-IR wavelengths. N 157B has 
been identified with SNR 0538-69.1, altho ugh its nature as a 
supernova remnant has been questioned llChu et al.ll2004l) . 
The age of the alleged s upernova remnant is only ~ 5000 yr 
( Wan g fc GottheUlll9 98) and it is therefore impossible to al- 
ready have led to the formation of stars in the nearby molec- 
ular cloud. It is possible that the formation of N 157BIRS1 
was triggered by feedback from the nearby OB association 
LH99, but we have no direct evidence for this. 

4.2 Triggered star formation in N 113 

N 113 is a smaller and perhaps simpler Hn region than the 
30 Dor mini-starburst complex. Fig. ll3l (bottom right) shows 
in great detail the interplay between the ionised gas and the 
molecular material in N 113. This image combines the Ha 
image (blue) with the 8/im IRAC/Spitzer image (red). We 
use here the pseudo-continuum image (Sect . |3~4*1 : this image 
contains not only stars (generally blue) but also most of the 
Ha emission structure; when combined with the IR image 
the ionised and neutral g well as the stellar content 

become apparent. This composite image presents a striking 
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example of the ionised bubbles created within the part of a 
molecular cloud facing massive early-type stars (Sect . UPl . 

In Fig. 1131 one can also see that the sequence of bright 
Ha knots A, B, E and G is co -spatial with the brightest 8/im 
emission. IWong et alJ i2006T) have recently performed aper- 
ture synthesis imaging of this dense region in N 113. The 
1.3-cm continuum observations reveal a string of 6 bright 
continuum sources, also located in the dense central region 
and a large molecular clump a few arcseconds away. Fig. 1131 
(bottom left) shows the 8/im Spitzer image of N113, with 
several sources identified. The positions of two water maser 
sources (lLazendic et al . 2002) are represented by circles. The 
east-most water maser source is the strongest detected in the 
Magellanic Clouds a nd it seems to be also asso ciated with 
an OH maser source jBrooks fc Whiteoaklll997l) and a rela- 
tively weak continuum source llWon^^t^Ll200al . The other 
water maser sour ce is associa t ed wit h the brightest contin- 
uum source that IWong et alJ i2006f) identify as a compact 
H n region. The approximate position of the dense molec- 
ular clump is also indicated in Fig. 1131 IWong et al] (120061) 
propose that the compact Hn region might trigger further 
collapse in this dense clump. 

The average molecular hydrogen density in the molec- 
ular cloud associated with N 113 is estimated at n c i ou d ~ 
200 cm" 3 llWong et al.l WoOG) which, for a cloud radius of 
Ncioud ~ 15 pc, corresponds to a column density of iVcioud ~ 
9 x 10 21 cm -2 . This just exceeds the threshold for a dif- 
fuse cloud to cool and collap se to for m stars, JVcri t icai ~ 
10 21 cm" 2 teergin et all20 04). However. IWong et al] 1120061 
argue that the cloud must be strongly clumped, which would 
imply that most (geometrically) of the cloud is much less 
dense — possibly iVcioud < iVcritical- The dense molecular 
clump that was detected within this cloud has a radius 
of -Rciump ~ 1.6 pc and a molecular hydrogen density of 
w c iu m p ~ 10 5 cm" 3 jWong et al.ll2006Tl : this yields a column 
density of N c i ump ~ 5 x 10 23 cm~ 2 , which comfortably ex- 
ceeds the threshold for star formation. The strongest den- 
sity enhancements are located near the interface with the 
hot gas; thus it is likely that the gas has been compressed 
as a result of feedback from the nearby OB stars. 

To summarise, the scenario that emerges in N113 is 
that in the compressed dense lane of neutral gas and dust 
star formation is occurring, pinpointed by the maser sources 
in its earliest stages and by the 1.3-cm continuum emission 
as the massive stars evolve. N 113 therefore presents a clear 
example of star formation triggered by the winds from mas- 
sive stars. A simila r scenario seems a lso likely for N 160A 
fFigJTSI see also lNakaiima et al]l2005|) . 




Figure 15. 8 fim IRAC/Spitzcr images of N 105A (top) and 
N 160A (bottom). In each im age, the scale indicates 1 arcmin. Cir- 
cles r epresent maser sources ISinclair et al]H992t 
2002), triangles are prot ostars (this work; 



Lazendic et al 



c prot( 

1984 1 ll^rmm^ et a .1 ll998tV and the square is a 



Epchtcin ct al 



. ... _ is a molecular 

core iBolatto et al] |2000D . A cross indicates the position of 
N105Ablob (Sect. l^2~2l . In both images North is to the 
top and East to the left and the greyscale is linear be- 
tween -20 and 250MJysr" 1 . The images are centered at ~ 
05 h 09 m 53 s , -68°53'27" and at ~ 05 ,i 39 m 41 s , -69°38'41" re- 
spectively for N 105A and N 160A. 



4.3 N 105A, a less evolved H 11 region? 

The situation seems different for N 105A. Fig. 1151 (top) 
shows the 8 /im IRAC/Spitzer image of N 105A, domi- 
nated by hot d ust and PAH emissio n, with the positions 
of N 105A I RS1 jEpchtein et al.lll984l: this work), water and 
OH maser iLazendic et al jl200a iBrooks fc WhiteoaJdll997b 
and methanol maser (ISinclair et al] I1992T) indicated. Al- 
though projection effects may affect our view, the mor- 
phology of N 105A shows little evidence for massive star 
feedback; no shell struct ure is seen in [O ill] images 
iAmbrocio-Cruz et aljH998T) , in spite of the proximity of the 
OB association LH31, identified near the molecular cloud 



iDopita et al 1ll994l) . Furthermore, there see ms to be no link 
betw e en the water and OH ma ser source dLazendic et al] 
120021 IBrooks fc Whiteoakl Il997l) and the protostar loca- 
tion and any e xternal trigger. How ever, the methanol maser 
to the south jSinclair et al.lll992)) does appear to lie in a 
dense knot at the rim of a cavity in the molecular cloud. A 
massive young stellar object, N105AIRS1, is already ion- 
ising the molecular cloud and also shows evidence of out- 
flows. It seems likely that it is only a matter of time before 
N 105 A IRS 1 starts sculpting the N 105 complex, influenc- 
ing current star formation occurring in the cloud core as 
signposted by the maser sources. 

The molecular cloud associated with N 105A is a factor 
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five denser than that in N113 llChin et aljll997l) . N 105A 
is therefore more likeiy to coliapse and fragment without 
the need for an external trigger. Thus N 105A seems to be 
at an earlier stage of its evolution, whereas N113, N 160A 
and in particular 30 Dor have seen more generations of stars 
forming prior to the current epoch of star formation. 



5 SUMMARY AND CONCLUSIONS 

We have conducted a survey for water maser emission 
in 30Doradus and a sample of Hn regions in the Large 
Magellanic Cloud, to investigate the conditions under 
which current star formation occurs. The locations of the 
maser sources are compared with infrared images from the 
ESO/VLT (at 3.8 pirn) and the Spitzer Space Telescope, and 
an Ha image in the case of N 113. We also present 3— 4/zm 
spectroscopy to investigate the nature of two protostar can- 
didates, in the N 157B and N 105A regions, showing evidence 
for the onset of the stellar feedback process in one of them. 

Our water maser survey of 30 Dor uncovered one new 
source, 0539— 691C, with a velocity consistent with the sys- 
temic velocity of the LMC. All detected water masers are 
located in the densest part of the nebula at the interface be- 
tween neutral and ionised gas. N 105AIRS1 shows strong H 
recombination line emission and its SED and IR colours are 
consistent with an embedded young massive star ionising its 
immediate surroundings. It also shows evidence for outflows. 
A nearby diffuse IR source, N 105A blob, shows both H re- 
combination and PAH emission. We identify N 157B IRS1 as 
an embedded protostar based on the analysis of its SED and 
a tentative detection of the 3.1 /im water ice feature. 

In the well-developed Hn regions 30 Dor, N 113 and 
N 160A, no water masers have been detected deep within 
the molecular cloud complexes. They are always found at 
the interfaces between molecular cloud and Hn region. This 
provides strong evidence that feedback from massive stars 
triggers subsequent star formation. Although in the dense 
cloud N 105A star formation seems to occur without evi- 
dence of massive star feedback, the general conditions in 
the LMC seem favourable for sequential star formation as a 
result of feedback. 

The wind speed and mass-loss rate of O and B stars 
decrease at lower metallicity (Mokiem & de Koter, in prepa- 
ration) , thus massive star feedback could be weak in metal- 
poor environments but we see no evidence for this in the 
LMC (metallicity ~40 per cent solar). This is important 
as the duration and intensity of star formation epochs in 
galaxies may depend on the efficiency of the local massive 
star feedback to trigger further star formation. 
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APPENDIX A: REGIONS WITH 
NON-DETECTIONS AT 22 GHZ 

We have observed several other Hn regions at 22 GHz, 
both in the Large and Small Magellanic Clouds. The po- 
siti ons listed in T ablc lA"T| had been previous l y obs erved 
by IWhiteoak et alJ (Il983l) and IScalise fc Brazl <1982Il We 
did not detect maser emission at any of these positions , 
confirming prev i ous n on-detections iWhiteoak et aljll983f) . 
lLazendic et al] i2002tl detected maser emission in N 159, 
with a peak flux density of 3.7 Jy at 247.5 k ms" 1 . No source 
was d etected in our spectrum. In the SMC, IScalise fc Brazl 
(1982) detected water maser emission in S 7 and S9 at the 
systemic velocity of the cloud (120 kms -1 ), with peak in- 
tensities of 7.4 and 4.2 Jy, respectively. We did not detect 
any water maser emission at either of these positions. Wa- 
ter maser emission is known to be intrinsically variable in 
gal actic star forming regions, both in in tensity and veloc- 
ity (iPalaei et al.lll993l : iTbfani et al.lll995f) therefore it is not 
unexpected that our survey does not detect all previously 
identified maser sources. 
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APPENDIX B: RDORADUS AND WHYDRAE 

The nearby galactic red giants RDoradus and WHydrae 
were observed at 22 GHz for reference. R Dor was observed 
several times each night. The resulting time series of the dif- 
ference with respect to the average spectrum is displayed in 
Fig.lBTl The average profile shows two distinct components, 
but only the high-frequency component shows variability. 
RDor pulsates radially in a semi-regular fashion with a pe- 
riod of P ~ 338 days, but the variability of the water maser 
emission occurs on a much shorter timescale — 9 days cor- 
responds to less than 0.03 cycles. 

The average spectrum is compared in Fie. lB2l to that 
observ ed 1.33 cycles previously (April 2000. Ivan Loon et alJ 
l2001bl) . The 200 1 spectrum looks more similar to the discov- 
ery spectrum ijLepine. Paes de Barros fc Gammon! Il97tl) . 
and is more symmetrical around the systemic velocity of the 
star as derived fr om the centroid of th e CO emission profile, 
v* = 24 kms -1 ijOlofsson et al J 1200 A . The double-peaked 
water maser profile suggests substantial radial amplification; 
this yields an estimate for the outflow velocity in the inner 
part of the dust envelope of «h 2 o ~ 2 kms -1 (peaks) to 
4 kms -1 (total extent), which is less than the outflow ve- 
locity in the outer part of the dust envelope as measured 
in CO, vco ~ 6 kms -1 . This implies that the wind is still 
accelerating in the region of the water maser. 

The water maser emission profile of WHya fFig. lB3j l 
has a single peak near the systemic velocity of the star, 
v* — 39.7 kms -1 (IGonzalez- Alfonso et alJll998T) . suggest- 
ing mainly tangential amplification. However, the emission 
profile has a broader pedestal component suggesting ra- 
dial motions of the order of vu 2 o ~ 2 kms -1 . The water 
maser is highly variable, but the maximum extent of the 
emis sion is well constrained and sugg ests vh 2 o < 4 km 
s" 1 dRudnitskii. Lekht fc Berulisl Il999)h As in the case of 
RDor, this is slower than the outflow velocity in the outer 
part of the wind as measured in CO, vco ~ 6.5 kms -1 
ilOlofsson et""al"ll2002t . and consistent with a wind that is 
accelerated through radiation pressure on dust grains. 

This paper has been typeset from a TpX/ I^TpX file prepared 
by the author. 
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Figure Bl. Time-series of the 22 GHz observations of RDor. The 
top panel shows each individual spectrum subtracted by the aver- 
age spectrum. The y-axis represents the actual time gaps between 
spectra. The average and variance spectra are shown in the panel 
below. Variability on the maser profiles seems to be restricted to 
the high frequency component at ~ 22 kms -1 . 
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Figure B2. Average 22 GHz spectra of RDor, from this work 
(full line) and Ivan Loon et al] i2001bl dot-dashed line). In 1.33 
pulsation cycles (see text), the maximum peak intensity has 
shifted from the low to the high frequency component. 
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Figure B3. 22 GHz spectrum of WHya. 



